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We have studied the transport properties of a large graphene double quantum dot under the 
influence of background disorder potential and magnetic field. At low temperatures, the evolution 
of the charge-stability diagram as a function of R-field is investigated up to 10 Tesla. Our re¬ 
sults indicate that the charging energy of quantum dot is reduced, and hence the size of the dot 
increases, at high magnetic field. We provide an explanation of our results using a tight-binding 
model, which describes the charge redistribution in a disordered graphene quantum dot via the 
formation of Landau levels and edge states. Our model suggests that the tunnel barriers separating 
different electron/hole puddles in a dot become transparent at high B-fields, resulting in the charge 
delocalization and reduced charging energy observed experimentally. 


I. INTRODUCTION 

Confining charge carriers in graphene continues to gen¬ 
erate interest owing to its customizable electronic prop¬ 
erties and compatibility with existing semiconductor de¬ 
vice processing PQ. Carbon atoms 12 C have low atomic 
weight and no nuclear spin (except for the 13 C isotope), 
so electronic interactions such as spin-orbit and hyper- 
fine coupling are expected to be weak in graphene, lead¬ 
ing to long electron-spin relaxation times [2j. Over the 
past decade lithographically-defined graphene quantum 
dots (GQDs) have proved to be an useful platform in 
which single electrons can be confined and manipulated. 
A number of experimental advances have been reported, 
such as charge detection [3], charge relaxation -3J, and 
electron-hole crossover [5] in graphene single quantum 
dots (GSQDs); and excited states [Bj, tunable interdot 
coupling :7j and charge pumping [8| in graphene dou¬ 
ble quantum dots (GDQDs). More recently, graphene 
quantum dots on hexagonal boron nitride (hBN) have 
enabled the influence of potential and edge disorder to be 
studied separately mm- Magnetic fields are a powerful 
tool for unveiling the nature of confined Dirac Fermion 
in GQDs. For example, the Fock-Darwin spectrum in 
the few-electron regime [5] and many-electron regime m 
as well as the Zeeman splitting of spin states m in 
graphene single quantum dot have been studied. On the 
other hand, while it is well-known that electron transport 
through graphene nanostructures is strongly affected by 
electron/hole puddles induced by potential fluctuations 
w\mm, detailed experimental and theoretical stud¬ 
ies are lacking to address this issue in GQD transport. 
In this letter, we study the effect of disorder by inves¬ 
tigating the transport properties of a large GDQD de¬ 
vice at magnetic fields in which Landau levels (LLs) are 
expected to form. At high enough B-field, our results 


suggest that electron/hole puddles in the dot tend to 
merge together, giving rise to a charge redistribution 
which can be observed experimentally. Our results are 
supported by tight-binding quantum simulations, which 
can be used to describe the charge redistribution in a 
disordered graphene quantum dot at high magnetic fields 
and gives deeper insight into our experimental data. 

II. COULOMB BLOCKADE MEASUREMENT ON 
A GRAPHENE DOUBLE DOT AT B=0 

Double quantum dots are a model system for inves¬ 
tigating the dynamics of electrons in a wide range of 
semiconductors mm- Charge stability diagrams - ob¬ 
tained by measuring the conductance as a function of the 
carrier density on each quantum dot - reveal a wealth 
of information about charging energy, interdot coupling 
and cross gate coupling strength, making them an ideal 
way to probe charge rearrangments in quantum dots at 
high magnetic fields. An atomic force microscopy im¬ 
age of the double quantum dot measured in this work is 
shown in Fig. [IJa). Our device consists of two litho¬ 
graphically etched (O 2 plasma) graphene islands each 
with a size of 200x250 nrn 2 labeled QDi and QD 2 in 
Fig. [IJa). They are mutually connected to each other by 
a 90 nm wide constriction and separately connected to 
the source/drain leads via two 80 nm wide constrictions, 
which act as tunnel barriers. Two plunger gates PG1(2) 
are used to tune the energy levels in QD^) while three 
side gates (SGI, SG2 and SG3) are used to tune the tun¬ 
nel barriers. The doped-silicon back gate (BG) is used 
to adjust the overall Fermi level. 

The measurements were performed in a dilution refrig¬ 
erator with an electron temperature around 100 mK. In 
Fig. [ljb), we show the measured differential conductance 
through DQD as a function of BG voltage (Vac= 20 /TV) 
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highlighting a region of suppressed current (the so-called 
transport gap [22] ) separating the hole from the elec¬ 
tron transport regime. At a back gate voltage within the 
transport gap (Vbg = 8.61 V, see arrow in Fig. 0 b )) we 
measure the DC current through DQD as a function of 
Vpax and Vpg 2 for a series of applied DC biases as shown 
in Fig. 0c) for 14=400 /TV, (d) for 14=1 mV and (e) for 
14=2 mV respectively. 

As expected, the current in the stability diagram 
evolves from triple points into bias-dependent trian¬ 
gles when the bias is increased. The horizontal and 
vertical measure of the honeycomb cell AVpgi and 
A Vp G 2 (Fig. Qd)) give the capacitances between the 
gate PG1 and QDi C s i~e/AVpGi=2.77 aF, and be¬ 
tween the gate PG2 and QD 2 C' S 2 ~e/AVpG 2 =l aF. 
Also the charging energies -Egi=0!i-AVpgi= 2.57 meV 
and Ec2 = ct-2-AVpG2=4-77 meV are obtained using the 
voltage-energy conversion factor CKi(2)=el4/<5VpGi(2) ex¬ 
tracted from the bias triangle as shown in Fig. Be). The 
difference in charging energies reflects the fact that the 
sizes of the dots are not equal and can be justified if the 
tunnel barriers defined by local disorder potential modify 
the size of GQDs mm- Within this picture, electrons 
from the source reservoir enter through a large localized 
state in QDi (Eci is small) to a small localized state 
in QD 2 (E C 2 is large), and then exit through the drain 
reservoir. Finally, the interdot coupling energy can be 
determined from the splitting of the triangles as shown 
in Fig. Be): E Cm =a 1 -AVp i G1 =0.29 meV. 


III. CHARGE STABILITY DIAGRAM IN 
PERPENDICULAR MAGNETIC FIELD 

The charge distribution in the QDs can be investigated 
by looking at how the charge stability diagram evolves 
under the influence of magnetic field. Fig. [2](a) shows the 
evolution of a region of the stability diagram, measured 
at T=100 mK, Vbg=8.61 V and I4=-l mV, for perpen¬ 
dicular magnetic fields ranging from 4 to 10 T. We have 
studied the stability diagram in a wide energy range and 
here only focus on four typical triple points for simplic¬ 
ity. The first observation is the field-dependent change 
in the dimensions of the honeycomb, which is highlighted 
by the dotted hexagonal outlines in Fig. [2ja) and is most 
pronounced from B =7 T to B =10 T, indicating the vari¬ 
ation in the capacitances C g \ and C g2 . In addition, the 
close-ups of the triangle in the bottom of each honey¬ 
comb, as shown Fig. [2j(b) , display a change in SVpci( 2 ) 
and implying that the conversion factors ai( 2 ) 

and interdot coupling energy Ecm also change with 13- 
field. It is worth noting that the size of triangle varies 
in the same honeycomb (i.e., B =7 T and 5=8 T), indi¬ 
cating the precise size of the localized state can change 
over a small range of gate voltage. The extracted charg¬ 
ing energy of each dot (Eci and Ec 2 ) and the interdot 
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FIG. 1: (a) Atomic force micrograph of the double quantum 
dot device measured in this work, (b) Measurement of the dif¬ 
ferential conductance through the DQD for varied back gate 
voltages. Data collected at V4c=200 /zV and T=1.4 K. Cur¬ 
rent through DQD as a function of Vpgi and Vpg 2 measured 
in a dilution fridge at T=100 mK with applied DC bias (c) 
Vsd= 400 pV, (d) Vsd= 1 mV and (e) Vso=2 mV. The po¬ 
sition of the triple points can be determined at low bias (see 
(c)), while at high bias the triple points evolve into triangles 
(see (d),(e)). 


coupling energy ( Ec m ) are shown in Fig. B c ) and (d) 
respectively. The QDs charging energies remain roughly 
unvaried (5 gi~ 3 meV, 5 g 2~5.5 meV) from 5=4 T to 
5=6 T then show a decreasing tendency from 5=7 T 
to higher fields. From 5=4 to 5=10 T, the percentage 
change in Ec 2 (37.3%) is larger than that in Eci (25%). 
By contrast, the interdot coupling energy shows a mono¬ 
tonic increase with field from 4 T to 10 T. Our results 
suggest that both dots increase their size at high 5-fields, 
which reflects on the decreasing charging energies. 


IV. MODEL AND SIMULATION 

It is well-known that the presence of charged impuri¬ 
ties in the Si0 2 substrate [231 12H or surface ripples [25] 
can induce electron-hole puddles with a size of tens of 
nanometres in exfoliated graphene flakes. This aspect 
considerably affects the electronic and transport proper¬ 
ties of graphene around the Dirac point and, as we will 
show, it plays a key role in our case. To take this into 
account, we consider a varying background potential V 
in a model QD, as shown in Fig. i a ) , where V fluctu- 
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FIG. 2: (a) The evolution of the charge stability diagram (taken at T=100 mK, Vbg=8.61 V and Vsd=- 1 mV) under the 
influence of perpendicular magnetic field from 4 to 10 Tesla, (b) Close-up of the triangle in the lowest position of each panel, 
as highlighted by the dashed square in the leftest panel of (a), (c) The charging energies of the QDs as a function of S-field. 
(d) The interdot coupling energy as a function of 5-held. 


ates from positive (red) to negative (blue) passing from 
V =0 (green). If V varies slowly, in each region of a large 
dot the energy bands will approximately correspond to 
the shifted energy bands of 2D graphene, as represented 
in Fig. |3](b) for zero (left panel) and high (right panel) 
magnetic field. At B =0 T, a gap is introduced to include 
the quantum confinement effects due to the dot. This 
gap progressively reduces at high 5-field along with the 
formation of Landau levels. Depending on the back-gate 
and background potentials, the Fermi energy E-p (here 
set to 0) and Dirac point can have locally different rel¬ 
ative positions, as indicated by dashed lines in the left 
panel of Fig. |3](b). The sign and strength of V deter¬ 
mine the nature (electron or hole) of the puddles and 
their density of state (DOS). We first consider the case 
5=0. For the F<<0 (F>>0) regions, the Fermi energy 
corresponds to level 1 (5) in Fig. [3jb) . As the level is far 
above (below) the charge neutrality point, it gives rise to 
the electron (hole) puddles with high DOS as shown in 
Fig- § c). In the region where F<0 (F>0), the Fermi 
energy corresponds to level 2 (4) and results in electron 


(hole) puddles with low DOS. In the region around V=0 
corresponding to the energy gap, the DOS is very low or 
0. These regions (the green region in Fig. [3jc) ) separate 
the puddles and can make the transport diffusive [26j . 

In the presence of high magnetic fields, due to the for¬ 
mation of LLs, part of the levels around the gap tends to 
the 0-th Landau level, thus reducing the gap. The other 
part rises, thus approaching the higher LLs as shown in 
the right panel of Fig. §b), where also the dispersive 
magnetic edge states are represented. In the low field 
regime, the LLs are far from being fully established and 
the DOS in the dot is low. However, the edge channels 
start developing with opposite chirality for electron and 
hole puddles as indicated by arrows in Fig. [3](d). At high 
magnetic field, the LLo is well developed with the con¬ 
sequent closing of the band gap. Therefore, in the F=0 
region the DOS is expected to increase and result in de¬ 
velopment of non-chiral channels connecting the puddles, 
as shown in the yellow region in Fig. |3j(e). At the same 
time, the other LLs start developing together with the 
chiral magnetic edge channels. In this regime, the DOS 
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FIG. 3: (a) Example of potential distribution in a large disor¬ 
dered quantum dot. (b) Left panel: Schematic band structure 
of a GQD in zero magnetic field. Right panel: in high mag¬ 
netic held. The red lines denote the electron-like levels and 
blue lines denote the hole-like levels. The black solid lines in¬ 
dicate the Dirac cone-likc dispersion, (c, d, e) Expected DOS 
distribution in the dot at zero held, low magnetic held and 
high magnetic held respectively. 


decreases in the bulk of the puddles while it increases at 
their edges. Electron transport through the dot is not 
confined in a particular puddle but can be delocalized 
in the dot through flowing in both the chiral edge chan¬ 
nels (red or blue arrows) and non-chiral channels (yellow 
region). 

In order to validate this picture, we performed numer¬ 
ical simulations of a QD with radius R= 47.5 nrn and a 
background potential consisting of two regions with V >0 
(81 meV) and V <0 (-66 meV), which determine the pres¬ 
ence of a hole puddle and an electron puddle, as shown in 
Fig. |4|a). The dot is described by a first-neighbour tight- 
binding Hamiltonian with a single p z orbital per atom 
and coupling parameter -2.7 eV. For the more details on 
the Green’s function formalism adopted for the simula¬ 
tions, refer to Ref. [27 . The calculated DOS of the dot 
is shown in Fig. |4](b)-(f) in arbitrary units. As expected, 
at low B (0 and 0.8 T) the DOS is low where F=0 and 
higher for larger \ V\. As B increases (2.8 T), the DOS de¬ 
creases a little in the centre of the electron/hole puddles, 
and it increases along the edge due to the progressive 
developing of magnetic edge states (Fig. id))- Note 
the presence of very high DOS region at the border of 
the dot. They correspond to zigzag edge sections, where 
very localized states appear 28]. At higher R-field, we 
observe the presence of high DOS in the V=0 region, 
which corresponds to the LL 0 , and the rise of edge states 
around the dot. The higher the field is, the larger the 
DOS is in the V=0 and edge regions, as can be seen for 
B>4.4 T in Fig. ie) and (f). 

The simulated background spectral current distribu¬ 
tion (which corresponds to the spatial distribution of 





FIG. 4: (a) Potential distribution in a quantum dot with 
R= 47.5 nm. The black line in the potential profile indicates 
the region where the potential is F=0. (b)-(f) Calculated 
local DOS in the dot at different magnetic field, (g)-(j) Cal¬ 
culated current distribution in the dot at different magnetic 
field. 


the conductive channels) [27j in the dot is shown in Fig. 
|4](g)-(j). At low magnetic field (5<0.8 T) the current is 
mainly concentrated in the high | V| regions and the V=0 
region seems to act as a barrier between the two puddles. 
At slightly higher field (B= 2.8 T) the current starts tend¬ 
ing to the V=0 region due to the progressive closing of 
the energy gap. At high fields (R>4.4 T), we observe 
the current flowing along the chiral magnetic edge states 
of the dot and along the non-chiral V=0 region, where 
the gap is now closed. In this regime, the current is de¬ 
localized in the dot and a charge rearrangemet can be 
seen compared to the case at low R-fields. Note that the 
more fractured the disorder potential is (meaning more 
existing electron/hole puddles) the more pronounced the 
charge delocalization effect will be at high fields. 
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V. DISCUSSION 

As the back gate voltage (Fbg= 8-61 V) where all the 
measurements were carried out is near the charge neutral¬ 
ity point, it is expected the background potential fluctua¬ 
tions will play a role and give rise to electron-hole puddles 
formed in the QDs [24] [29]. In this situation, our model 
can be readily adapted to explain the data. Due to the 
closing of the energy gap at high enough magnetic fields, 
at the F=0 region the DOS is high and develops non- 
chiral channels which connect the puddles, as sketched 
in Fig. [3je) . Hence, the current can flow through the 
puddles via crossing the high DOS non-chiral channels 
at the interface, thus making electrons no longer con¬ 
fined in a particular puddle but delocalized in a larger 
puddle, resulting in a smaller charging energy. Here we 
point out that the field has to be high enough for the LLs 
and edge states to be fully developed to close the energy 
gap. The threshold 5-field for this to happen in a GSQD 
with a relatively smaller size (50 nm diameter) is around 
10 T [5]. As a result, the change in charging energies in 
our case is most pronounced from 5=6 T to 5=10 T for 
QDi and from 5=7 T to 5=10 T for QD 2 , as shown in 
Fig-d c). The threshold 5-field for QD 2 is higher owing 
to its larger charging energy (smaller puddle), in which 
the magnetic length f's=y / ?i/eB has to be comparable or 
even smaller than the puddle size. The magnetic length 
for 5=7 T is around 9 nm, implying the size of puddle in 
the dot is around (or more than) twice the critical mag¬ 
netic length, in good agreement with the puddle size (20 
nm) measured in graphene B51 . In addition, we observed 
the change in Eq 2 from 5=4 T to 5= 10 T is larger than 
that in Eci ■ This is expected, since Eqi is smaller than 
Ec 2 , indicating electrons tunnel through a larger local¬ 
ized state in QDj and a smaller localized state in QD 2 . In 
other words, the disorder potential in QD 2 is more frac¬ 
tured than that in QDi. Therefore, at high 5-field the 
charge delocalization effect is more pronounced in QD 2 
than that in QDi, giving rise to the larger 5-dependent 
charging energy in QD 2 . Here we note that the GQD has 
to be large for the substrate disorder to play an important 
role, which may be the reason that the decreasing charg¬ 
ing energy with 5-field is not observed in other relatively 
smaller GQD with large edge-to-bulk ratio Bn or GQD 
on hBN pTO] where substrate disorder is less important. 
The increasing interdot coupling energy may be also un¬ 
derstood as the charges rearranging from the center of 
the puddles to the edge of the dot (see Fig. @(g)-(j))- 
This scenario depends on the progressively formed edge 
state with increasing 5-field, and can be observed in the 
whole range of 5-field (Fig. [2jd)) . 

A more recent study of GNR on hBN substrate has 
indicated that the localized states may also extend into 
the leads of the device, giving rise to smaller charging 
energies than expected from the geometry of GNR alone 


[501 . However, two conditions are crucial for this effect 
to be seen. One is the substrate disorder has to be much 
weaker than the edge disorder, and the other is the edge- 
to-bulk ratio of device has to be large enough for the 
edge to play an important role. Transport in GQDs on 
hBN is dominated by the edge roughness for QDs with 
diameter less than 100 nm [2], Each condition is met by 
their relatively small GNR (30 nm x 30 nm) on hBN, 
thus the edge disorder is strong enough to localize elec¬ 
tron wavefunction along the edge to the leads. On the 
contrary, our large dot (and the tunnel barrier GNR) 
with smaller edge-to-bulk ratio should diminish the in¬ 
fluence of the rough edges on overall transport, meaning 
localization along the edge still happens but transport 
is dominated by bulk contributions. Therefore, we ar¬ 
gue that the charge redistribution (based on substrate 
disorder) in our QDs is the main factor that leads to a 
variation of dots area and contribute to the decreasing 
charging energies in magnetic fields. 

The effect of disorder can be also seen in Fig. [5j where 
a stability diagram measured in different cool-down of 
device at 5-fields (a) 3.2 T, (b) 3.8 T and (c) 4.4 T is 
presented. The triangle shape first distorts at 5=3.2 T 
and then splits into two separated ones (Fig. [ 5 ] (b)), then 
moves further apart and form an additional row of tri¬ 
angles (Fig. [ 5 ] (c)). We attribute this newly appeared 
triangles to the formation of a localized state in a mag¬ 
netic field, which is capacitively coupled to the original 
dots gU02]- A schematic illustration is shown in Fig. [ 5 ] 
(d) and (e) to address such a scenario. When a localized 
state is formed in magnetic field, while the gate voltage 
is sweeped it can add or subtract charges discretely to 
the parasitic dot, thus altering the entire environment 
abruptly and unexpectedly. The fact that the splitting 
occurs on both gate spaces suggests the localized state 
can affect two dots in a similar way, implying its location 
is in the central GNR (Fig. [ 5 ] (d)). The new dot acts as a 
gate which will shift the triple-points in charge stability 
diagram, consequently, leading to an additional row of 
triangles adding adjacent to the original ones. 

VI. SUMMARY 

In summary, we have fabricated and studied the 
magneto-transport properties of a large GDQD device. 
In different cool downs, we observed a honeycomb pat¬ 
tern which is typical of charge stability diagrams for a 
DQD system. We studied the evolution of the charge 
stability diagram under the influence of 5-field up to 10 
T. The charging energy and the interdot coupling energy 
show different dependence with 5-field, suggesting the 
size of both dots become larger in high field. Our inter¬ 
pretation is supported by numerical simulations in which 
we show the confined charges in the puddles of GQDs can 
be redistributed from the bulk to the edge through the 
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FIG. 5: The charge stability diagram measured at Vbg= 9 V 
and 14=1 mV in (a) B= 3.2 T, (b) B= 3.8 T and (c) 13=4.4 T 
showing a formation of an additional dot under the influence 
of magnetic field and strong couple to the original dots, (d) 
Graphic illustration of an effect of a localized state formed in 
magnetic fields, (e) Same as (d) but with a charge added into 
the localized state. It capacitively couples to the original dots 
and forces the DQD to reconstruct its wavefunction. 

formation of LLs and edge states. At high enough 13-field, 
due to the closing of energy gap, electrons are delocalized 
via crossing the non-chiral channels connecting different 
puddles, resulting in a smaller charging energy. 

This work was financially supported by the EPSRC. 
Additional data related to this publication is available at 
the xxxxxxx data repository. 
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